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ABSTRACT 

I compare the density profile of dark matter (DM) halos in cold dark matter (CDM) 
N-body simulations with 1 h _1 Mpc, 32 h^ 1 Mpc, 256 h _1 Mpc and f 024 h" 1 Mpc box 
sizes. I compare the profiles when the most massive halos are composed of about 10 5 
DM particles. The DM density profiles of the halos in the f h" 1 Mpc box at redshift 
z ~ 1 show systematically shallower cores with respect to the corresponding halos in 
the 32 h _1 Mpc simulation at z ~ 3 that have masses, Md m , typical of the Milky Way 
and are fitted by a NFW profile. The DM density profiles of the halos in the 256 h _1 
Mpc box at z ~ are consistent with having steeper cores than the corresponding 
halos in the 32 h _1 Mpc simulation, but higher mass resolution simulations are needed 
to strengthen this result. Combined, these results suggest that the density profile of 
DM halos is not universal, presenting shallower cores in dwarf galaxies and steeper 
cores in clusters. More work is needed to validate this finding at z = 0. Physically 
the result sustains the hypothesis that the mass function of the accreting satellites 
determines the inner slope of the DM profile. But the result can also be interpreted 
as a trend with the dynamical state in the assembly process of halos of different mass. 
In comoving coordinates, r, the profile 
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Pdm oc 



X a (l + x) 3 - a 

with X = c A r/r A and a = (9 + 3n)/(5 + n) w 1.3 + {M^ u - 1)/(M$£ M + 1), 
provides a good fit to all the DM halos from dwarf galaxies to clusters at any redshift 
with the same concentration parameter ca ~ 7. Here, r A (Md m ) is the virial radius, n 
is the effective spectral index of the initial power spectrum of density perturbations 
and Mdm,u — Mdm/(3 X 10 14 M Q ). The slope, 7, of the outer parts of the halo 
appears to depend on the acceleration of the universe: when the scale parameter is 
a = (1 + z)^ 1 ;$ 1, the slope is 7 w 3 as in the NFW profile, but 7 w 4 at a > 1 
when £1a ~ 1 and the universe is inflating. The shape of the DM profiles presents a 
significant scatter around the mean. It is therefore important to analyse a significant 
statistical sample of halos in order to determine the mean profile. 

I compare the DM profiles in the 1 h" 1 Mpc box with the same simulation in- 
cluding stars, baryons and radiative transfer presented by Ricotti, Gnedin and Shull. 
Radiative feedback effects produce a larger scatter in the density profile shapes but, 
on average, do not affect the shape of the DM profiles significantly. 

Key words: cosmology: theory, dark matter - galaxies: dwarf, clusters, formation, 
halos - methods: numerical, N-body simulations 



1 INTRODUCTION 

According to the currently favoured galaxy formation sce- 
narios, galaxies are formed from the gravitational growth 
of tiny density perturbations imprinted on a uniform uni- 
verse during inflation. Small mass perturbations grow faster 
and constitute the building blocks for the assembly of larger 
galaxies and clusters. In the nonlinear phase of the gravita- 



tional collapse the dark matter (DM) particles are shock 
heated to the virial temperature and settle into a dark 
halo with mean overdensity about A = 178 (according to 
the simple spherical collapse model in a fiat firj = 1 uni- 
verse). Analysi ng N-body simulations of hierarchica l struc- 
ture formation iNavarro. Frenk fc White) (Il996tll997ft (here- 
after, NFW) have proposed that DM halos develop a univer- 
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sal density profile, valid for virialized halos at any redshift 
and with any mass, from dwarf galaxies to clusters. Fur- 
ther works have shown that an approximatively universal 
profile develops regardless of the details of the adopted cos- 
mology and initial power spectru m of density perturbations 
<Huss et alll999UEke et alJ200ll) . The NFW density profile 
has a core cusp p<xr _1 and in the outer regions the density 
decreases as p oc r~ 3 . The details of the gravitational col- 
lapse affect only the relative size of the core with respect to 
the virial radius (the halo concentration parameter, ca) but 
not the slope of the profile. This result implies that during 
the virialization process the memory of the mass function 
and profile shapes of the building blocks that formed the 
halo is lost. 

Recent observations of dwarf galaxies (spheroidal and 
irregular) and low surface brightness (LSB) galaxies seem 
to show that a flat core for the DM density profile 
is more compa tible with the measureme n ts of the ro- 
tatio n curves (|v^n_dgji^c«ch_^t_alj [gjjOOj; |de_Fjtok_^t_alJ 
| 200ll: Ide Blok fc Bosmal 120021: iBorriello fc Saluccil l200lt 
ISalucci fc Burkertll200[F . These observations have renewed 
interest in the subject and a number of solutions have been 
proposed to solve this possible discrepancy. The proposed so- 
lutions range from feedback effects of the baryons and stars 
on the DM density profile to modifications of the physi- 
cal properties (e.g., scattering cross section, temperature) 
of the DM particles. S ome examples of a lternatives to CDM 
are warm dark matter feode et all200ll) and self-interacting 
DM JSpereel fc SteinhardtfeOOoTT The potential discrepancy 
between the steep DM cores found in N-body simulations 
with the observed flat core of dwarf galaxies is perhaps the 
most serious problem faced in CDM cosmologies today. 

High resolution simulations and careful convergence 
tests have been done or are underway to u nderstand if 
the discrepa ncies found by some authors fe.g.. lMoore et alJ 
( 1998. fl999T) find steeper core profiles) are due to numerical 
artifacts. In the core of DM halos the particle crossing time 
is much shorter than the Hubble time and the integration 
of trajectories is subject to su btle numerical effec ts that are 
difficult to keep under control. IPower et alJ (I2002T) estimated 
that about 1 million particles are needed in order to have 
convergent results in the inner 1% of the halo virial radius. 
To have such a high number of particles per halo the use 
of a special technique is required. This technique consists of 
simulating with high mass resolution only the particles that 
will end up in the halo of interest at z = 0, while having less 
mass resolution for the particles that end up far away from 
the halo of interest. This has the drawback that each simu- 
lation can resolve only one halo at a time, therefore selection 
effect biases and a poor statistical sample could affect the 
reliability of the final result even if the simulation is very 
accurate and has high resolution. Note that most of the de- 
tailed and computationally expensive work on the subject 
has focused on simulating halos at z — with typical mass 
similar to the Milky Way. 

The mass function of the building blocks of dwarf galax- 
ies is different in comparison to the one for larger galaxies 
or clusters. On small scales (or masses) the mass fraction 
of virialized halos is about constant as a function of the 
logarithmic DM mass of the halos because the power spec- 
trum has a slope n ~ — 3. This means that small and large 
mass satellites contribute equally to the mass of the accret- 
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Figure 1. Linear power spectrum of density fluctuations at z = 
(solid line) and logarithmic slope, n(k), of the power spectrum 
(dashed line). The horizontal bars show the range of wave num- 
bers, k = 2n/\, resolved in the 1,32,256 and 1024 h" 1 Mpc 
simulations. 

ing halo. In massive halos as in the Milky Way or in clusters 
of galaxies, the contribution from large mass accreting satel- 
lites is instead dominant. For this reason it is more crucial 
to resolve very small mass satellites (i.e., to have high mass 
resolution) in simulating dwarf galaxies than in simulating 
massive galaxies or clusters. 

In this paper I compare the density profiles of DM ha- 
los in four CDM N-body simulations that are identical apart 
from the box sizes that are 1 h _1 Mpc, 32 h _1 Mpc, 256 h _1 
Mpc and 1024 h _1 Mpc. The simulation method is a two- 
level particle-particle particle-mesh (P 3 M) with 256 3 parti- 
cles and ACDM cosmology. In dimensionless units the only 
difference between the simulations is the initial power spec- 
trum of density perturbations. I compare the profiles when 
the most massive halos are composed of about 10 5 DM par- 
ticles and the clustering in the simulations is similar (see 
Fig. |5J . This is done in order to minimise possible system- 
atic errors in the simulations and to be more confident that 
any variation in the profile shapes is produced by the differ- 
ent initial power spectrum. The masses of the halos studied 
in the four simulations, from the small to the large box, are 
typical of dwarf galaxies, Milky Way mass galaxies, clusters 
and superclusters of galaxies, respectively. 

If the cores of small mass halos are dense enough to 
resist disruption and survive "undigested" when incorpo- 
rated into a larger object, they should determine the halo 
structure in the inner regions. Us ing scaling arguments 
ISubramanian. Cen fc Ostrikerl i2000T) have shown that, in a 
flat universe with an initial power spectrum of density per- 
turbation P(k) oc k n , the core density profile is p& m oc r a 
with a = (9 + 3n)/(5 + n). The numerical results that they 
report approximately confirm this expectation. In ACDM 
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cosmologies the logarithmic slope of the initial power spec- 
trum ranges between —3 < n < 1, depending on the mass 
scale (see Fig.0. From dwarfs to superclusters the logarith- 
mic slope of the core density profile should range between 
^ a 2. As shown in Fig.0 the power spectrum in the 1 
h _1 Mpc box is a power law with n ~ —2.7. According to 
the aforementioned scaling relationship a ps 0.4. This slope 
is in agreement with what I find in the 1 h -1 Mpc box. In the 
32 h^ 1 Mpc box the power spectrum is not a single power 
law but the slope is close to n ~ —2 and therefore a ~ 1 
is expected, consistent with the NFW universal profile and 
with the result found for this simulation. In cluster and su- 
percluster size halos a should be larger than in galaxy size 
halos (a ~ 1.5). Indeed, in the 256 h -1 Mpc simulation I 
find slopes of the inner profile a ~ 1.4. But this is not a 
strong result because the number of particles in each halo 
is barely sufficient to start seeing a deviation with respect 
to the NFW profile. Moreover, the introduction of the pa- 
rameter a for the slope of the inner part of the halo profile 
(which can be easily calculated knowing the mass of the 
galaxy halo) does not increase the number of adjustable pa- 
rameters in the profile fitting formula. This happens because 
the concentration parameter, C/[ FW (Mdm), needed to fit the 
profiles with the NFW formula, is in this case a constant. 
The changing slope of the inner profile mimics the variation 
of the concentration parameter in the NFW profile, shifting 
the value of the radius where the circular velocity reaches 
its maximum value (see appendix [XJ. 

This paper is organised as follows: in § [5] I explain the 
methods for the simulations and in §|3]I show the results for 
dwarf galaxies, normal galaxies and clusters of galaxies. In 
§2] I demonstrate how a density profile with changing inner 
slope a can fit all the profiles from dwarfs to clusters with a 
constant concentration parameter. The discussion and con- 
clusions are presented in §|S| The equations for the circular 
velocity and integrated mass for the NFW profile and for 
a halo with arbitrary inner and outer slopes of the density 
profile are contained in appendix 1X1 



NFW universal profile. The differences of the profiles in this 
two simulations are likely to be caused by the different initial 
power spectrum of density perturbations, since the they have 
identical smoothing length in dimensionless units, identical 
Fourier modes of the initial density perturbations, similar 
number of integration time steps and comparable number 
of particl es per halo. According to this paradigm and the 
results of ISubramanian et alJ J2000l) . the slope of the DM 
profiles in simulations with box sizes larger than 32 h _1 Mpc 
should be steeper than the NFW profiles. I perform two ad- 
ditional simulations with 256 h _1 Mpc and 1024 h _1 Mpc 
box sizes to test this hypothesis. The large box simulations 
present the problem that in order to achieve the same degree 
of clustering as in the smaller boxes it is necessary to evolve 
the simulation in the future when a ^ 1 and Ha > 0.7. For 
the 256 h _1 Mpc simulation the required level of clustering 
is reached when 0.9 < a < 1.5 and the number of particles 
per halo is barely sufficient to start noticing a steepening 
of the inner profile with respect to the NFW profile. In the 
case of the 1024 h" 1 Mpc box the level of clustering reached 
in the smaller boxes can never be achieved because the uni- 
verse begins to inflate (the scale factor becomes a > 1000 in 
a few time steps) and the event horizon starts decreasing, 
preventing any further clustering. I present the results for 
this simulation because they are interesting for understand- 
ing what determines the outer slope of the density profile. 
The number of particles in these halos is too small for mea- 
suring the slope of the inner profile but it is evident that, in 
an accelerating universe, t he outer slope of the densit y pro- 
file is steeper. Analytically. ISubramanian et alJ i2000T) show 
that the slope of the outer profile is 7 ~ 4 in a low density 
universe and in an accelerating ACDM universe. Table 
lists the scale factors when the simulations with different 
box sizes, Lbox, have run for a comparable number of time 
steps and the most massive halos are composed of the same 
number of DM particles. 

In the following paragraphs I explain in greater detail 
the numerical techniques and the sanity checks adopted to 
derive the density profiles and the rotation curves. 



2 THE METHOD 

In the literature most efforts have focused on studying 
galaxies at z = w ith masses typical of the Milky Way. 
iRicotti et all (l2002allbh have performed high-resolution sim- 
ulations of the formation of the first galaxies using a cos- 
mological code that solves the equations for the DM parti- 
cles, baryons, stars and radiative transfer. These simulations 
achieve a mass resolution of M p = 4.9 x 10 3 Mq for the DM 
using 256 3 particles and 1 h _1 Mpc box size. In these simu- 
lations at z ~ 10 or, expressing the redshift in terms of the 
scale factor, at a = (l + z)" 1 ~ 0.09, the profiles of the most 
massive halos (Mdm, ~2x 10 s Mq) appear flatter than the 
NFW profile. I therefore perform the same simulation for 
only the DM particles, finding again systematically flatter 
halo cores than found by NFW. Since this result could be 
produced by subtle numerical problems of the N-body inte- 
gration or the small number of particles in the halo (typically 
between 10 4 and 10 5 ), I repeat the simulation for a 32 h _1 
Mpc box down to a = 0.3 when the most massive halos have 
the same number of DM particles as in the 1 h' 1 Mpc box 
at a = 0.09. In this case the profiles are well fitted by the 



2.1 N-body Simulations 

I use an N-body code based on the P 3 M method with 
two levels of mesh refinement feertschinger fc Gelblll99ll : 
iGnedin fc Bertschingej |l996). In all the simulations the 
number of particles is 256 3 and the Plummer softening pa- 
rameter is e = 0.02 in cell units AL, where AL = Lbox/256 
and Lbox is the box size. Unless stated otherwise, I use 
scale-free (dimensionless) units in this paper. The cosmology 
adopted is a flat ACDM model with Q A = 0.7, Q, m = 0.3 
and h — 0.7. The initial conditions for the density and 
velocity fields are calculated using the COSMIC package 
jBertschingerill9 95") assuming a scale-invariant (tilt n — 1) 
initial spectrum of density perturbations and ACDM trans- 
fer function. The spectrum, shown in Fig. is normalised 
imposing a variance erg = 0.91 in spheres with radius of 8 
h _1 Mpc at z = 0. The initial conditions are calculated us- 
ing the Zel'dovich approximation until z = 166,91,48 and 
26 for the L bo x = 1,32,256 and 1024 h" 1 Mpc boxes, re- 
spectively. The N-body simulations start at those redshifts. 
I use the same random realization for the wave phases and 
directions in all simulations. 
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Table 1. Clustering and time steps in the four simulations, a is the scale factor; N p has two columns: the numbers on the left show the 
number of particles in the most massive halo, the right contains the mean number of particles for the ten most massive halos. 
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The halos a re identified using DENMAX 

feertschinger fc Gelbl Il99ll) with smoothing parameter 
G = 1000. The DENMAX algorithm identifies halos as 
maxima of the smoothed density field with smoothing 
length Lq — Lbox/G. The algorithm assigns each particle 
in the simulation to a group (halo) by moving each particle 
along the gradient of the density field until it reaches 
a local maximum. Unbound particles are then removed 
from the group. The results of DENMAX depend on the 
degree of smoothing used to define the density field. A finer 
resolution in the density field will split large groups into 
smaller subunits and vice versa. This arbitrariness in results 
is a common problem of any group-finding algorithm, and 
it is not only a numerical problem but often a real physical 
ambiguity. Especially at high redshift, since the merger rate 
is high, it is difficult to i dentify or define a si ngle galaxy 
halo. In appendix B in iRicotti et all (|2002eJ| the mass 
functions obtained using DENMAX with different smooth- 
ing lengths are compared with the analytical expectation 
using the Press-Schechter formalism. It was found that the 
smoothing parameter G = 1000 gives the best results for 
boxes with 256 3 particles. 

In each simulation with a different box size the pro- 
files of the halos are compared when the level of clustering 
is similar and the positions of the most massive halos co- 
incide. This happens after a comparable number of time 
steps, but the number of steps required becomes larger as 
the box size increases since the halos are more concentrated. 
The most massive halos in the simulations are composed of 
about N p = 10 4 — 10 s DM particles and the typical number 
of time steps to reach this level of clustering is a few thou- 
sands (see Table 0. Fig. [5] shows the sizes and positions of 
bound halos in the 1 h _1 Mpc box at a = 0.085 and in the 
32 h" 1 Mpc box at a — 0.230, when the clustering in the 
two simulations is similar. 



2.2 Density Profiles 

The density profiles are calculated in three different ways: 

(i) Using the list of particles identified as halo members by 
DENMAX. This means all the bound particles. 

(ii) Using, additionally, the particles identified by DENMAX 
as unbound particles. 

(iii) Using all the particles. This includes all the halo satel- 
lites and neighbouring halos. 

Usually profiles using (i) and (ii) are almost indistinguish- 
able. But the profiles using (iii) look different in the outer 
parts if the halo is disturbed by accreting satellites. The pro- 
files obtained with method (iii) are used for the analysis as 
in previous studies. The other profiles are useful to discrim- 
inate between isolated halos and halos that are undergoing 
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Figure 2. Position of DM halos in the 1 h _1 Mpc box at a = 
0.085 (left panels) and 32 h -1 Mpc box at a = 0.23 (right panels). 
The circles show the sizes (virial radii) and positions of the 100 
most massive halos in the simulation. The 10 most massive halos 
are shown as thicker circles. The positions of the most massive 
10000 halos are shown as points. 



accretion of massive satellites. One can choose whether to 
exclude or include in the analysis halos that are undergo- 
ing major mergers; either choice does not change the result 
significantly. 

The profiles are shown in dimensionless units: the radii 
are expressed in number of cells x (the cell size is AL = 
Lbox/256 in comoving coordinates) and the masses in num- 
ber of DM particles N p . The mass and density of each halo 
are sampled as a function of the radius, x. in shells with 
constant logarithmic spacing dx/x = 0.01. The dimension- 
less circular velocity is more appropriate to compare the 
profiles since the virial radii are difficult to determine when 
the halos are not isolated. Moreover, the relative error on 
the density profile is smaller at radii corresponding to the 
maximum of the circular velocity since the number of parti- 
cles per shell at intermediate radii is maximum. Therefore, 
if the halo profile is not disturbed by accreting satellites, 
the location and the value of the maximum circular veloc- 
ity provides the best method to compare halos of different 
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Table 2. Scaling constants. 
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a 
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3.9 


3.69 


7.26 


0.07-0.09 


32 
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8.21 


235.7 


0.23-0.3 


256 


1000 


10.9 


1885 


0.9-1.5 


1024 


4000 


12.7 


7542 
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masses. The dimensionless circular velocity is defined as, 

_ V c (x) _ / N dm (x) 
Vc ~ ~VT ~ V ~x~ 

where V c ° = (GM P /AL) 1/2 = 7.365 x \Q- 2 (L box /l ft^Mpc) 
km s _1 . If the slope of the inner density profile is p(r) oc r~^, 
the circular velocity has a slope V c oc r' 1- *'' 72 - 1 . The Pois- 
son error for the density is Sp/p = dN p (x)^ 1 ^ 2 , where 
dNp(x) is the number of particles in each shell of radius 
x. The error on the mass as a function of the radius is 
SMdm/Mdm(x) = N p (x)~ 1 ^ 2 , where N p (x) is the number 
of particles in each sphere of radius x. The error on the cir- 
cular velocity is 5v c /v c = 0.5SMdm/Mdm. Table [5] lists the 
values of useful constants for converting dimensionless units 
into physical units. 



3 RESULTS 

In Figs. 13151 (top) the dimensionless circular velocities are 
shown for the nine most massive halos in the simulations 
with box sizes L box = 1, 32, and 256 h~ Mpc, respectively. 
The statistical errors, Sv c , on each point are about the 
sizes of the symbols and therefore are not shown to avoid 
excessive crowding. I show the profiles obtained using 
methods (ii) (small dots) and (iii) (large dots). I do not 
show those obtained using method (i) since they are almost 
indistinguishable from (ii). I show the fit with a NFW profile 
obtained by matching the position and value of the maxi- 
mum of v c and the virial radius (arrow) for the halos that 
are not severely perturbed by mergers. The points are con- 
nected in the range of radii that are reliable: x m i < x < xa- 
Here x A = 2.37(A/178) _1/3 (iVp/10 4 ) 1/3 is the dimension- 
less virial radius and x m % is the radius that contains a 
number of particles N p (x m i) 300(p(x m i)/10 4 ) 1//2 . This 
condition ensures that two-body rel axation is not impor tant 
in flattening the inner profile (see iPower et al.l (I2002T) for 
details on convergence studies). Note that x m i is always 
larger than twice the Plummer softening length 2e = 0.04, 
therefore an incorrect force calculation does not affect 
the results at those radii. In summary, the sanity checks 
adopted are as follows. 

1) Only radii well above twice the softening length are 
considered. 

2) Only radii that enclose N p ^ 300(p/10 4 ) 1,/2 particles are 
considered. 

3) No merging halos: profiles (ii) and (iii) give similar- 
results for the value and radius of the maximum circular- 
velocity. 

4) The same criteria are adopted for all the simulations 
with different box sizes, Lbox, an d scale factor a. 

Figs. I3TB1 (bottom) show the mean profiles obtained by 



scaling (shifting in logarithmic scale) the radii to the same 
Xmax = 1 and normalising the maximum of v c to v c nax = 1. 
Each radial bin, logarithmically spaced, is weighted equally 
so long as N p ^ 300(p/10 4 ) 1/2 . Only the halos that are 
not undergoing a major merger are considered. If all the 
ten profiles or a different subset of profiles are considered, I 
find that the resulting mean profile does not change signif- 
icantly. Also, using the arithmetic or geometric mean does 
not change the mean profiles. A possible criticism to this 
method is that the halos analysed are not isolated as in 
works where a single galaxy is re-simulated with higher res- 
olution. Time-dependent effects could thus produce flatter 
or steeper cores. For this reason, I present the mean profile at 
different redshifts (listed in the figure captions). If the halos 
are not relaxed to their final configuration the mean profiles 
should differ at different redshifts. Moreover, the increased 
statistical sample of halos makes the result more robust. 

In the 1 h _1 Mpc simulation (see Fig. |2J the slope of 
the inner profiles is, on average, flatter (i.e., the circular ve- 
locity rises more steeply) than the NFW profile. The typical 
masses of the halos in this simulation are Mdm ~ 10 s Mq. 
Even at radii close to v max the circular velocity starts to 
deviate from the NFW fit. The slope of the outer profile 
is instead consistent with 7 = 3 as in the NFW formula. 
Note that method (ii) does not attribute several DM parti- 
cles to the halo, producing an outer profile that decreases 
more steeply. It is not clear if DENMAX fails not to in- 
clude those particles or if they effectively do not belong to 
the halo. The definition of halo boundary (or virial radius) 
is somewhat arbitrary and time dependent since halos are 
constantly growing. Especially for small mass halos the flat 
mass function of the accreting satellites and the short ac- 
cretion timescale makes it more difficult to define the outer 
edge of the galaxy. This is also true for dwarf galaxies at 
2 = 0, even if the accreting satellites are mostly invisible 
because they do not contain stars. For this reason an "iso- 
lated" dwarf galaxy cannot be found or defined in the same 
way as for massive galaxies for which the mass function of 
accreting satellites is dominated by halos with masses sim- 
ilar to the host halo. There are halos for which the use of 
either method (iii) or (ii) gives a different radius where the 
circular velocity is maximum. Those halos are excluded from 
the analysis since they are severely disturbed by accreting 
satellites (but if they were to be included in the analysis, 
the mean slope would remain flatter than a — 1). The mean 
density profiles are best fitted by inner slopes in the range 
q ~ 0.4 — 0.5, equivalent to mean slopes of the inner circular- 
velocity (3 ~ 0.75 - 0.8. 

In the 32 h" 1 Mpc simulation (see Fig.QJ the NFW pro- 
file provides a good fit to all the halos that are not severely 
disturbed by accreting satellites. The typical masses of the 
halos in this simulation are M dm ~ 5 x 10 12 M Q . The mean 
density profiles are best fitted assuming inner slopes in the 
range a ~ 0.9 — 1, equivalent to mean slopes of the inner 
circular velocity f3 ~ 0.5 — 0.55. 

In the 256 h _1 Mpc simulation (see Fig. |KJ the inner 
density profile slopes appear slightly steeper than predicted 
by NFW (i.e., flatter circular velocities). Unfortunately, the 
number of particles in the inner region is barely sufficient 
to start seeing a discrepancy. A larger number of particles 
per halo is needed to construct reliable rotation curves at 
small radii where the discrepancy, expected by extrapolat- 
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Figure 3. (Top). Dimcnsionless circular velocity, ti Cl as a func- 
tion of comoving radius, x, for the nine most massive halos in the 
1 h _1 Mpc simulation at a = 0.085. The points are connected in 
the interval where the profile is reliable (the arrow indicates the 
virial radius). The smaller points are obtained using method (ii) 
(see text) and the thin solid line shows the NFW fit to the profile 
(for those profiles that are not undergoing major mergers). (Bot- 
tom). Normalised mean circular velocity v c (top panel) and its 
slope (bottom panel) as a function of the normalised radius x of 
the nine more massive halos in the simulation. The thick dashed 
line shows the NFW profile and the thick solid line a profile with 
a = 0.2. The mean profiles at at a = 0.070, 0.080, 0.085 and 0.090 
are shown to ensure that the result is not affected by time de- 
pendent effects such as mergers and to increase the statistical 
significance of the result. 
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Figure 4. (Top) Same as in Fig. EJfor the 32 h" 1 Mpc box at 
a = 0.23. (Bottom) Same as in Fig.|3]for the 32 h" 1 Mpc box at 
a = 0.23, 0.25, 0.27 and 0.30. 



ing the reliable part of the rotation curve, should become 
evident. The typical masses of the halos in this simulation 
arc M dm ~2x 10 15 M Q . The mean density profiles are best 
fitted by inner slopes in the range a ~ 1.3 — 1.4, equivalent 
to mean slopes of the inner circular velocity f3 ~ 0.3 — 0.35. 
In Fig.\5\(bottom) , it appears that the slope of the outer pro- 
file becomes steeper when a > 1. For this reason, included 
in the figure are the best fits for two cases: (i) when the 
outer slope of the profile is 7 w 3 at a = 1.0, 1.2, and (ii) 
when the slope is 7 m 4.0 at a — 1.5, 1.7. This dependence 
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Figure 5. (Top) Same as in Fig.EJfor the 256 h" 1 Mpc simula- 
tion at a = 1.0. (Bottom) Same as in Fig.[3]for the 256 h _1 Mpc 
box at a = 1,1.2,1.5 and 1.7. The solid and dotted thick lines 
show two profiles both with a = 1.4, but with slopes of the outer 
profile 7 = 3 and 7 = 4, respectively. 



of the outer profile slope on the scale factor becomes more 
evident in the 1024 h^ 1 Mpc simulation and is discussed in 
the following paragraph. 

In Fig. |S] (top), the dimensionless circular velocities are 
shown for the nine most massive halos in the simulation with 
box size L^ ox — 1024 h _1 Mpc at a = 10. The number of 
particles in these halos is too small for measuring the slope 
of the inner profile but it is evident that the outer slope 
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Figure 6. (Top) Same as in Fig.EJfor the 1024 h" 1 Mpc simu- 
lation at a = 10. (Bottom) Density profile (solid lines) and slope 
of the density profile (dashed lines) for the 1024 h _1 Mpc sim- 
ulation at a = 10. The thin solid line is the fit with the NFW 
profile. The outer slope of the density profile is 7 « 4, steeper 
than in the NFW profile. This is probably caused by the acceler- 
ating universe. 



is steeper than the NFW profile. Fig. |S| (bottom) shows the 
density profiles (solid lines) and slope of the density profiles 
(dashed lines) for the nine more massive halos in the same 
simulation. The outer density profile has a slope 7 ~ 4. At 
a — 2.3 the slope is 7 ~ 3.5 and increases to 7 ~ 4 when 
the scale factor is a — 4, 10, 100 and 1000. The physical 
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explanation of this result is likely to be related to the ex- 
pansion rate of the universe. When Qa{z) ~ at a $ 0.5 
the universe is decelerating and 7 ~ 3. As a increases to 
a ^ 1 the universe starts inflating (at a = 2, £Ia ~ 0.95 
in a flat universe) and the slope of the outer profiles be- 
comes 7 ~ 4. It is possible that even at 2 = the outer 
profile of recently virialized clusters is slightly steeper than 
7 = 3. According to the aforementioned interpretation, the 
observation of a steep outer profile in nearby clusters would 
in dicate that the universe is a ccelerating (see also the results 
of ISubramanian et alJ feOOOT) for low density universes). In 
§ 0] I show that for the halo profiles in this simulation the 
radii, r max , where the circular velocity is maximum, decrease 
with increasing scale factor, consistent with the predictions 
assuming an inner slope of the density profile a ~ 1.5. 

3.1 Radiative Feedback Effects on the DM profile 

Fig.[7|shows the mean v c for the nine most massive halos in 
the simulation 256Llp3 presented bv lRicotti et alJ (|2002b) 
at a = 0.075 and a = 0.090. This simulation attempts to 
simulate realistically the formation of the first galaxies in 
the universe modeling radiative feedback effects in the early 
universe. The parameters of the simulation are the same as 
in the 1 h _1 Mpc simulation presented in this work. But 
in addition to DM particles it includes gas dynamics, star 
formation using a Schmidt-Law, metal enrichment from star 
formation and radiative transfer. Radiative feedback effects 
produced by UV radiation emitted by the first stars are cal- 
culated following the molecular processes involving H2 for- 
mation, destruction and cooling. From the results presented 
in this work it appears that radiative feedback processes are 
not responsible for the flattening of the DM profile in small- 
halos. The slopes of the inner profile of the halos in this 
simulation are similar to the corresponding N-body simula- 
tion, or perhaps slightly steeper. 



4 CONSTANT CONCENTRATION DM 
PROFILES 

In comoving coordinates, r, all of the simulated halos in this 
work can be fitted by a profile with the form (c/., Zhao 1996) 

pdm{x) x X^l + X)^ ' (1) 

where X — CAr/rA, )*a is the comoving virial radius, a = 
ct(Md m ) is a function of the halo mass, Mdm, and ca ~ 7 is 
the concentration parameter that, as I will show, is a con- 
stant. In physical coordinates, R, it is simply X — caR/Ra, 
where Ra — (1 + z)~ 1 r& is the virial radius. Given the ini- 
tial power spectrum of density perturbations, the parameter 
q is a function of the halo mass only. The values found in 
this w ork agree with the predictions of ISubramanian et alJ 
( 2000) (see Fig.|HJ|. An acceptable fit to the theoretical value 
of a for the case of the ACDM cosmology, adopted in this 
work, has the form a » 1.3 + (M^ 14 - l)/(M^f >M + 1), 
where Mdm,i4 = Md m /(3 x 10 14 M©). This fit is five percent 
accurate for masses 10 6 < M dm /Mq < 10 18 . 

The logarithmic slope, 7, of the outer parts of the halo 
appears to depend on the acceleration of the universe: when 




Log x 

Figure 7. Same as in Fig. |3] (Bottom) for the 1 h _1 Mpc simu- 
lation presented by Ricotti et al. (2002b) at a = 0.075 and 0.090 
(see text). 

the scale parameter is a ,$ 1, the slope is 7 « 3 as in the 
NFW profile, but is 7 « 4 at a > 1 when Qa ~ 1 and the 
universe is inflating. Equation Q is normalised imposing the 
condition that the integrated mass inside the virial radius 
is m d m{X = ca) = M dm (see appendix |A"1. The comoving 
virial radius is defined as 

,A.(32kpc)fA) fJ^V 73 , 

to satisfy the relation Mdm ~ (47r/3)poArA, where A is the 
mean overdensity of the halo and po is the mean DM density 
at 2 = 0. A depends on the cosmology since it is related to 
the mean overdensity and the time it takes for a density 
perturbation to stop expanding with the universe and "turn 
around," starting the collapse. According to the simple top- 
hat spherical collapse approximation, A « 1 78^0(2)° ,45 for a 
flat universe with a cosmological constant l|Eke et alJ [l998L 
I now explain why this profile has the property of a con- 
stant concentration parameter. The changing slopes of the 
inner and outer parts of the profile mimics the variation of 
the concentration parameter in the NFW profile, shifting the 
value of the radius, r m ax, where the circular velocity reaches 
its maximum value. For the NFW profile (i.e., for a — 1) 
TA/rmax = Ca FW /2.16, therefore the location of r ma x with 
respect to the virial radius is inversely proportional to the 
concentration parameter. In the general case 

= CA ( 2 ) 
r ma x f(a, 7) 

The location of the maximum circular velocity depends on 
the inner and outer profile slopes and on the concentration 
parameter. In appendix^ I show that a good approxima- 
tion for the function f(a, 7) is / ^ (A/j) 2MS (2 - a). In 
Fig. H it can be seen that the values of rA/r ma x (circles), 
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Figure 8. r^/r max (circles) and a (triangles) for the simulations 
listed in Table [I] as a function of the halo DM mass M^ m . The 
open circles show r^/r max when 7 > 3 at a > 1. The dashed line 
shows the theoretical value of a (see text). The solid and dash- 
dotted lines show r/^/r max for the profile given in Eq. using 
7 = 3 and 7 = 4, respectively. 

for all the simulations in this work (listed in Table 0, agree 
with equation @, assuming constant concentration param- 
eter ca = 7. Moreover the values of a (triangles) are in 
go od agreement with the th eoretical relationship proposed 
bv lSubramanian et alJ feOOd) . a = (9 + 3n)/(5 + n) (dashed 
line), where n is the effective slope of the power spectrum. 
The solid and open circles show r&/r m am measured in pro- 
files with 7 = 3 (i.e., at a < 1.2) and with 7 > 3, respec- 
tively. The solid and dash-dotted lines show r^jtmax given 
by equation using the theoretical value of a (shown by 
the dashed line), 7 = 3 and 7 = 3.5 — 4, respectively. In the 
1024 h" 1 Mpc box simulation I find that 7 ~ 4. The inner 
slope of the profile cannot be measured because of the small 
number of particles in the inner regions. But the measured 
values of rA/r mal are consistent with values of the inner 
slope a ~ 1.5, predicted by the theory. 



5 DISCUSSION 

The results presented in this work and their interpretation 
seem to disagree with most previously published work, al- 
though the simulations presented here are not of substan- 
tially higher resolution than those in the literature. In this 
section I analyse the discrepancies with previous works and 
point out that the results do not disagree strongly with most 
previous works, save for the 1 h _1 Mpc simulation. The re- 
sult of the 1 h _1 Mpc simulation cannot be compared with 
previous works since a similar simulation has not been per- 
formed before. The joint analysis of the three simulations 
with Ltox = 1, 32 and 256 h^ 1 Mpc suggests that the slope of 



the inner part of the DM profile depends on the power spec- 
trum of linear perturbation. This interpretation of the re- 
sults disagrees with the common view of universality of DM 
profiles, although several groups have found results that con- 
tras t with this idea, b oth in the context of g alactic dynam- 
1CS (iNipoti et alJl200l e.g. ,) and cosmolo gy (|SveTfc_Whitej 
Il99ct iKravtsov et aljll99a Ijing fc Sutoll2000l : I Jind 120001) . 
Here I identify and summarise three possible reasons why 
the present results differ from a large number of published 
work on N-body simulations. 

(i) The profiles of halos with typical masses resolved by the 
1 Mpc simulation have not been analysed before. This 
simulation is equivalent to a scale free simulation with linear 
power spectrum P(k) tx k n with n = —2.7 and Einstein-de 
Sitter cosmology, since at z > 2 the cosmological constant 
contribution is negligible. Previous works have studied scale 
free si mulations with power spectrum index n ^ — 2. For in- 
stance ICole fc Lacevl (11996) studied the n = 0, — 1, — 2 cases 
and lNavarrc^^dTlu99^) the n = 0,-0.5,-1,-1.5 cases. 
The simulations in those works have the same resolution as 
in the present study. Their results do not contrast with the 
findings of this work because when n ^ — 2 the inner slope 
of DM profiles is also a ~ 1 — 1.5. But, as mentioned, the 
resolution of the simulations (about 10 4 — 10 5 particles per 
halo) is marginally sufficient to discern between a ~ 1 and 
q ~ 1.5. The n = —2.7 case is markedly different from the 
others since it corresponds to a scenario in which all mass 
perturbations become non-linear at about the same time 
(an index n < —3 would correspond to the top-down for- 
mation scenario, in which small mass halos form after the 
large ones). For this reason, I believe that the result found 
analysing the 1 h" 1 Mpc simulation is new and complements 
previous studies. 

(ii) Interpretation of the results. Theoretically it is un- 
clear what determines the slope of the inner DM pro- 
files. Some works argue that the slo pe depends on the 
power spectrum of initial perturbations llSver fc Whitell99ci 
ISubramanian et alJl200ol) . others argue that t he slope tends 
asymptotically to a stee p slope a ~ 1 — 2 iAlvarez et alJ 
l2003t iDekel et al] 1200 j) . Recent still unpublished high- 
resolution N-body results from Navarro and collaborators 
show that the inner slope does not converge to 01 = 1 but 
the slope slowly flattens to a £ 1. This result contrasts with 
most aforementioned semi-analytical results indicating that 
physical processes that determine the density profile of DM 
halos and their time evolution are not well understood. 

The halos analysed in the present work are the most mas- 
sive halos at a given redshift. For this reason they are re- 
cently virialized halos with different masses, compared after 
about the same number dynamical times from their forma- 
tion. The halos analysed in high-resolution N-body studies 
are usually chosen to be "isolated": halos that did not ex- 
perience "recent" major mergers. Therefore, depending on 
their mass, their age from the virialization differs greatly at 
z — 0. Semi-analytical works study the equilibrium profile 
of halos given an accretion history. The different evolution- 
ary states of the halos studied in literature could explain 
the discrepancies in the results if the time scale to reach the 
equilibrium configuration is long compared to the accretion 
time. I find that the halo profile depends on the accretion 
history. But it is possible that this is true only during a tran- 
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sient period of time, before the halo relaxes to a different 
equilibrium configuration. This concern is especially valid 
for the 1 h" Mpc simulation since in a scale free simulation 
with n ~ —2.7 the halos continuously accrete satellites. The 
result is therefore statistically significant at redshift about 
10 but the shallow slope of the profiles might not survive 
to redshift z = 0. Even if the shallow profiles are transient, 
depending on the time scale to reach the equilibrium config- 
uration, the result at z = 10 might suggest a reason for why 
not all dwarf galaxy hal os at z = presen t shallow cores. 
According to the work of lDekel et al.l (I2003f) . substructure is 
not tidally striped but instead is compressed when the halo 
profile is shallow. Therefore substructure could survive undi- 
gested for quite a long time. Moreover even at z — 0, small 
mass halos, due to their larger number density, are expected 
to capture dark satellites of non-negligible mass with higher 
frequency than larger mass halos. Further work is nee ded, 
and is currently in progress (iRicotti fc Wilkinsonll2003) . to 
extend the result of the present work from redshift z = 10 
to redshift z — 0. The aim of this study is to understand 
if the flat profiles at z = 10 are in a stable configuration 
that would survive for a Hubble time and if the profiles are 
consistent with the observed flat cores of dSph galaxies in 
the local group. Preliminary results show that the dynami- 
cal state of the halos is stable and the inner profiles remain 
flat for 20 Grys if the halo is evolved in isolation. 

To further support the case for a dependence of the inner 
profile on the mass of the halos I show that flat inner profiles 
are associated with an almost isotropic velocity dispersion of 
DM particles while steeper profiles have a velocity dispersion 
slightly biased in the radial direction. This result is shown in 
Fig- El The lines show the parameters (3 = l — ag/a 2 (dashed 
lines) and /3* = 1 — (vg)/(v 2 ) (solid lines) as a function of the 
scaled radius x = r/r max , where r max is the radius where 
the circular velocity is maximum. Here, a'g = {v 2 } — (vg) 2 and 
a 2 — (v 2 ) — {v r } 2 , where v r is the radial velocity and vg is the 
tangential velocity. The /3 parameter measures the degree 
of anisotropy of the velocity dispersion: (3 — corresponds 
to an isotropic velocity dispersion, f3 < to a tangentially 
biased velocity dispersion and < /3 < 1 to a radially biased 
velocity dispersion. The lines from top to bottom show the 
mean value of (3 and /3* for the 9 most massive halos in the 
256 h" 1 Mpc, 32 h" 1 Mpc and 1 h" 1 Mpc simulations. The 
shallow inner DM profiles in the 1 h _1 Mpc simulation have 
almost isotropic velocity dispersion. Steeper DM profiles in 
the 32 h _1 Mpc and 256 h _1 Mpc simulations have a velocity 
dispersion slightly biased in the radial direction. This result 
is consistent with the analytical expectations based on the 
solution of the Jeans equations in spherical systems. 

(iii) Method and analysis of the results In this paper I 
present statistical results without preselecting the halos to 
re-simulate with higher resolution as has been done in pre- 
vious studies. In order to compare simulations with obser- 
vations a statistically significant sample should be analysed. 
Halos that did not reach an asymptotic equilibrium con- 
figuration (if they exist) or halos presenting substructures 
produced by recent mergers should be included in the anal- 
ysis if they appear to be statistically significant. Especially if 
there is not an easy way to determine the dynamical stage of 
a halo because the accreting substructure is invisible (small 
halos at z — should accrete mostly dark halos), under- 
standing what fraction of dark halos of a given mass is in 




Log x 



Figure 9. (3 = 1- <J 2 e lo 2 (dashed line) and (3* = 1 - (vj) / (v 2 ) 
(solid lines) as a function of the scaled radius x = r/r max . The 
lines from top to bottom show the mean value of (3 and (3* for 
the 9 most massive halos in the 256 h — 1 Mpc, 32 h — 1 Mpc and 
1 h _1 Mpc simulations (shown in the top panels of Figs. 3-5). 
The shallow inner DM profiles in the 1 h — 1 Mpc simulation have 
almost isotropic velocity dispersion (f3 ~ 0). Steeper DM profiles 
in the 32 h — 1 Mpc and 256 h — 1 Mpc simulations have a velocity 
dispersion slightly biased in the radial direction (i.e., ay > (Tg). 



a transient configuration is relevant to solve the problem 
of the DM cores observed in a fraction of galaxies at low 
redshift. 

As expected the halos studied in this work show substan- 
tial substructure. 3D images of the halos show that numer- 
ous satellites are present, independently of the halo mass. 
Perhaps in small mass halos substructure can survive "undi- 
gested" closer to the halo cores. Nevertheless, in most cases, 
this does not affect the approximately spherical symmetry of 
the halos and the determination of the halo centres. The halo 
centres are calculated, using DENMAX, as the minimum of 
the gravitational potential. I also calculate the centres as the 
maximum of halo density and as the centre of mass of bound 
particles. The halo centres calculated using those alternative 
definitions are practically coincident and choosing either of 
them does not change the profiles. Exceptions are profile 8 
in the 1 h _1 Mpc simulation and profile 4 in the 256 h _1 
Mpc simulation that have maximum density centres offset 
by 2e, where e is the Plummer softening parameter. 

In this work the inner slope of the profiles and the core 
radii are determined fitting the circular velocities instead 
of the density profiles. This procedure partially remove the 
degeneracy between the slope a of the inner profile and the 
core radius r 3 . This is illustrated in Fig. 1101 which shows 
the maximum likelihood estimate of the parameters a and 
r„ fitting the density profile (left) or the circular velocity 
(right). The la, 2a and 3a confidence contours are shown 
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assuming constant relative errors on the density and circular 
velocity in each logarithmically spaced radial bin. 



6 CONCLUSIONS 

Two main problems faced by CDM cosmologies are re- 
lated to the properties of small mass galaxies. Namely, (i) 
the number of visible galactic satellites (dwarf galaxies) 
in the Loca l Group is smaller than predicted by N-body 
simulations jMoore et alJll99St iKlvpin et alJ|l999l) and (ii) 
the flatness of the DM cores in dwarf and LSB galaxies 
is not reproduced by N-body simulations. The first prob- 
lem can be solved by including radia tive feedback effects 
JChiu et alJl200ll: iRicotti et aill2002bl) in cosmological sim- 
ulations. Galaxies with masses M dm £ 10 8 - 10 9 M are too 
small to retain the gas that is photoevaporated by stars or 
by the ionising background after reionization. Consequently, 
the luminosity of most dwarf galaxies is predicted to be too 
low to be detected (i.e., the mass-to- light ratio should in- 
crease for halos with smaller masses) or zero for very small 
mass halos. Here I have shown results of N-body simulations 
suggesting that the second problem might not contradict the 
predictions of CDM cosmologies. But further investigation is 
needed to understand if the dwarf galaxies showing shallow 
DM profiles at z ~ 10 (e.g., just after the time when most 
of the dwarf galaxies formed) are in a configuration that 
will survive unchanged until redshift z = 0. In summary the 
main findings of this work are: 

• The slope, a of the inner profile of DM halos is deter- 
mined by the mass function of the accreting substructure 
(i.e., the initial power spectrum of initial perturbations). 
Dwarf galaxies have on average flatter DM cores, and clus- 
ters, steeper cores than galaxies similar to the Milky Way 
for which the NFW profile is a good fit to the halo density 
profile. 

• The logarithmic slope, 7, of the outer parts of the halo 
appears to depend on the acceleration of the universe: when 
the scale parameter is a ;$ 1, 7 « 3 as in the NFW profile, 
but 7 ~ 4 at a > 1 when Qa ~ 1 and the universe is inflating. 

• A density profile in the form p d m tx X~ a (l + X) ,+a , 
where X — cat/ta, a is a function of the halo mass 
(a » 1.3 + [M^ 14 - 1]/[M^ 6 14 + 1], where M dmM = 
M dm /[3 x 10 14 M ]) and 7 ~ 3 (but 7 ~ 4 if fl A ~ 1), can 
fit all the profiles from dwarfs to superclusters with constant 
concentration parameter ca — 7. This provides a physical 
explanation for the variation of the concentration parame- 
ter in the NFW profile as a function of the halo mass: the 
changing slopes mimic the variation of the concentration 
parameter by shifting the radius where the circular velocity 
reaches its maximum value. 

• The dependence of a on the mass of the halo, and there- 
fore on the slope of the power spectrum, n, agrees with the 
theoretical expectation a = (9 + 3n)/(5 + n), derived assum- 
ing that "undigested " satellites determine the halo structure 
in the inner regions iSubramanian et alJl20o6T) . Note that if 
the DM halo develops a flat core, the chance f or accreting 
satell ites to survive tidal disruption is larger JPekel et alJ 
120031) . 

The shape of the profiles for a given mass and red- 
shift present significant statistical variations. This can be 



attributed to two factors. Firstly, galaxies are not isolated 
entities and especially at high redshift are severely disturbed 
by neighbouring satellites or by undergoing major mergers. 
This produces irregularities in the azimuthally averaged den- 
sity profile. Secondly, the accretion history of satellites has 
statistical variation and depends on the local environment. 
A dwarf galaxy forming at high redshift from a large a peak 
of the power spectrum of initial perturbations (the type of 
galaxy simulated in this work) could have an inner profile 
different from the same mass galaxy forming at lower red- 
shift from a smaller a perturbation. Also, the slope of the 
profile could differ if a galaxy forms in a void or in an over- 
dense region, even for galaxies with similar masses and radii. 

I have identifi ed three main reasons to explain w hy pre- 
vious studies fe.g.. lEke et alJ200lHMoore et al Jl999ll found 
that the inner slope of the DM profiles are not determined 
by the power spectr u m or accretion h istory (but see also 
ISver fc White] Jl998ft : iKravtsov etafl (Il99gl) : Ijing fc Sutol 
feOOOl) : IJind fcOOCf) for a different point of view), (i) The 
profiles of halos with typical masses resolved by the 1 h _1 
Mpc simulation have not been analysed before, (ii) The dy- 
namical state of the the halos at z — 10 might be different 
at z = 0. But there are not published high-resolution simu- 
lations for dwar f-size halos (Mdm ~ 1 8 Mr) at z = 0. The 
classic results of lCole fc Lacevl il996l) ; iNavarro et alJ (Il997f) 
are based on the same type of simulations presented here, 
(iii) Most of the computational efforts have focused on sim- 
ulating a single galaxy-size or cluster-size halo with mass 
Mdm ~ 10 12 Mq at z = with higher resolution (about 10 6 
particles per halo) than in this study. Such a large number 
of particles is required to study the profile in the very in- 
ner regions of the halo. But the task of achieving reliable 
profiles with such a high spatial resolution is sensitive to 
numerical inte gration erro rs and requires careful resolution 
studies llPower et alJl2002T) : perhaps this is a reason for the 
disagreement between groups, using different codes, on the 
slope of the inner profile. In this work I find that the shape of 
the DM profiles present significant scatter around the mean. 
It is, therefore, important to analyse a significant statistical 
sample of halos in order to measure the mean profile. Works 
where a single halo is re-simulated with higher resolution are 
affected by a bias, difficult to control, determined by the cri- 
teria for picking the halos to re-simulate. Finally, the fitting 
method could be important as well for the results. The de- 
generacy between the slope of the inner profile and the value 
of the core radius is more easily broken fitting the circular 
velocities instead of the density profiles. This is because the 
radius where the circular velocity is maximum is univocally 
determined by the value of the core radius, r s (for the NFW 
profile is r max = 2.16r s ). The halo profiles at any redshift 
and mass, when renormalised matching the values and radii 
of the maximum circular velocities, have to be identical if 
their profiles can be all fitted by the NFW profile. 

iMoore et alJ lll999l) find steep profiles in cluster size ha- 
los (a ~ 1.5), in agreement with the results of this work. The 
results of this work also agree with the c lassic NFW result for 
Milky Way size galaxies but not with lMoore et all <ll999l) for 
halos with the same mass. In this work the spatial resolution 
is not as large as in the aforementioned studies, therefore I 
cannot study the slope of the DM profile in the very centre 
of the halos. Nevertheless, I find evidence for disagreement 
with the NFW predictions even at relatively large radii. Per- 
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Figure 10. Maximum likelihood estimate of the core radius, r s , and inner slope, a. The la, 2a and 3a confidence contours are shown. 
The degeneracy between the core radius and inner slope when the density profile is fitted with the generalised profile given in Eq. 
(left panel) can be partially broken fitting the circular velocity (right panel). 



haps this new result has been found because it is the first 
time that a simulation of halos with masses typical of dwarf 
galaxies and mass resolution M p ~ 10 3 Mq has been per- 
formed. The result for the cluster size halos could also be 
understood with similar reasoning (i.e., the cluster masses 
are larger than in previous studies: Mdm ~ 10 15 — 10 16 M@). 
The steeper core density profile found in clusters may not 
be as evident as the flatter core found in dwarf galaxies but, 
combined with the result for normal galaxies, the findings 
from the three simulations strengthen the case for an inner 
slope that changes with the mass of the halos. 
I conclude this work with two final remarks: 

- According to the scaling relation a — (9 + 3n)/(5 + 
n), the inner slope of DM profile of dwarf galaxies is more 
sensitive to the power spectrum index, n, than more massive 
galaxies: Sa — 1.58n if n ~ —3. A tilt 8n ~ 0.2 in the initial 
power spectrum should produce a variation Sa « 0.3 of the 
inner slope of the DM profile. 

- The best place to look for the effects of the cosmological 
constant is the outer profile of clusters and superclusters at 
low redshift. The prediction is for a slope steeper than 7 = 3. 
Perhaps gravitational lensing studies could be useful in this 
regard. 
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APPENDIX A: GENERALISED PROFILE 

The comoving virial radius, r&, here is defined as the ra- 
dius where the halo mean overdensity is A, i.e., Mdm = 
(47r/3)poAr^. Here Mdm is the mass of the halo, po is the 
mean DM density at z = and A « 178 if Q(z) = 1. A 
density profile with inner slope a, and outer slope 7, can be 
written in the form 

, \ _ PoA£(c, a , 7 ) 
Pdm(r) - XQ(1+X)7 _ a 

where X — r/r s , r s — ta/ca is the core radius and <5( ClQ , l7 ) 
is a normalisation constant. The NFW profile has a = 1 and 
7 = 3. The mass enclosed inside the radius, X, is 

rX 



If 



dy 



(i + v)* 1 



V 



47rp A(5( c , clj7 )r 3 



9 



(Al) 



2Fi(g;'y - a;g + 1;-X) 



where g = (3 — a) and 2F1 is the hypergeometric function. 
If 7 = 3, the hypergeometric function simplifies as 



2 F 1 (g;g;g+l;-X)=gY, 



B(n;g) 
(9 + n)n 



(-A)", 



where B is the beta function. When a = 1 and 7 = 3 (NFW 
profile) equation lAH can be integrated analytically and the 
solution is 



I v\ 47T : NFW 3 

m dm (X) = —p A5 c r s 



ln(l + A) 



X 



l + X 
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to be evaluated numerically. The simplest fit to the numeri- 
cal solution, shown in Fig. lAlf b). is given by / = (4/j) 2 ' 68 . 
Using this fit in conjunction with the previous for f(a, 3), a 
general fitting formula is found to be 



/(«, 7 )«(4/ 7 f 68 (2-a), 



(A3) 



Figure Al. (a) f(a, 7) = r max /r B as a function of a for 7 = 3 
(circles), 7 = 3.5 (open circles), and 7 = 4 (squares). The solid 
lines show the fit to the numerical solution, (b) f(a = 0,7) = 
Tmax/Ts as a function of 7. The solid lines show the fit, / = 
(4/7) 2 65 (2 — a), to the numerical solution. 



The normalisation constant S 



(c,a ,7) 



is such that irid m (X 



ca) — M dm . For the NFW profile, 5£"' w = c A /[ln(l + c A ) - 
ca/(1 + ca)]- In the general case equation 1AH has to be 
evaluated numerically. If 7 = 3, analytical solutions exist for 
a = 0,0.5,1,1.5 and 2: 



3 

x c 
5 



[_(i+™ +ln(1 + x)], if Q = 

+ ArcSmh(A 1/2 )] , if 



3 



(1+4X/3)X 1/2 



ln(l + X) 

X 



X 
1+X 



+ ArcSmh(A 1/2 

if a = 2 

If 7 = 4 an analytical solution exists for a generic a 
1 + X 



if a = 1.5 



2 F 1 (-X) 



X 



[l-3(l+X)- 1 +3(l+X)- 2 -(l+X)- 



The circular velocity is defined as 

Vc{X) : 



m dm (X) 



Xr s 



(A2) 



For the NFW profile equation HA2H has a maximum V™ ax oc 
ta at r m ax = 2.16r s . For a fixed value of r A , V™ ax depends 
weakly on ca, and r max is inversely proportional to ca- In 
the general case, r max = /(ct, 7 )r s and must be calculated 
numerically. In Fig. lAlf a) I show the numerical solution 
and a simple fit to /(a, 7) = r max /r s as a function of a for 
7 = 3, 3.5 and 4. If a = the maximum, /, of v c can be 
found solving the equation 



2[l-(l + /) 7 ] + (7-l)/ 



(7 - iff + 7/ + 



27 



= 0. 



7-1. 

If 7 = 4 the solution is / = 2, but in general the solution has 



which is at least 5 percent accurate in the parameter interval 
< a < 2 and 3 < 7 < 4. 



REFERENCES 

Alvarez M. A., Aim K., Shapiro P. R., 2003, in M. Reyes 
E. V.-S., ed., "The Eighth Texas-Mexico Conference on 
Astrophysics" RevMexAA SC 

Bertschinger E., 1995, COSMIC, GC-3 report, (astro- 
ph/9506070) 

Bertschinger E., Gelb J. M., 1991, Computers in Physics, 
5, 164 

Bode P., Ostriker J. P., Turok N., 2001, ApJ, 556, 93 
Borriello A., Salucci P., 2001, MNRAS, 323, 285 
Chiu W. A., Gnedin N. Y., Ostriker J. P., 2001, ApJ, 563, 
21 

Cole S., Lacey C, 1996, MNRAS, 281, 716 

de Blok W. J. C, Bosma A., 2002, A&A, 385, 816 

de Blok W. J. C, McGaugh S. S., Bosma A., Rubin V. C, 

2001, ApJ, 552, L23 
Dekel A., Devor J., Hetzroni C, 2003, MNRAS, in press 

(astro-ph/0204452) 
Eke V. R., Navarro J. F., Frenk C. S., 1998, ApJ, 503, 569 
Eke V. R., Navarro J. F., Steinmetz M., 2001, ApJ, 554, 

114 

Gnedin N. Y., Bertschinger E., 1996, ApJ, 470, 115+ 
Huss A., Jain B., Steinmetz M., 1999, ApJ, 517, 64 
Jing Y. P., 2000, ApJ, 535, 30 
Jing Y. P., Suto Y., 2000, ApJ, 529, L69 
Klypin A., Kravtsov A. V., Valenzuela O., Prada F., 1999, 
0.5 ApJ, 522, 82 

Kravtsov A. V., Klypin A. A., Bullock J. S., Primack J. R., 

1998, ApJ, 502, 48 
Moore B., Governato F., Quinn T., Stadel J., Lake G., 

1998, ApJ, 499, L5 
Moore B., Quinn T., Governato F., Stadel J., Lake G., 

1999, MNRAS, 310, 1147 
Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 

563 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 
493 

Nipoti C, Londrillo P., Ciotti L., 2003, MNRAS, accepted 
Power C, Navarro J. F., Steinmetz M., 2002, The inner, 

(astro-ph/0201544) submitted 
Ricotti M., Gnedin N. Y., Shull J. M., 2002a, ApJ, 575, 33 
Ricotti M., Gnedin N. Y., Shull J. M., 2002b, ApJ, 575, 49 
Ricotti M., Wilkinson M. I., 2003, in preparation 
Salucci P., Burkert A., 2000, ApJ, 537, L9 
Spergel D. N., Steinhardt P. J., 2000, Physical Review Let- 
ters, 84, 3760 

Subramanian K., Cen R., Ostriker J. P., 2000, ApJ, 538, 
528 

Syer D., White S. D. M., 1998, MNRAS, 293, 337 
van den Bosch F. C, Robertson B. E., Dalcanton J. J., de 
Blok W. J. G., 2000, AJ, 119, 1579 



© 2002 RAS, MNRAS 000, EE1 



14 M. Ricotti 

Zhao H., 1996, MNRAS, 278, 488 



© 2002 RAS, MNRAS 000.H1IT41 



